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Abstract 

In order to gain insight into the possible Ground State of Quantized Einstein's 
Gravity, we have derived a variational calculation of the energy of the quantum 
gravitational field in an open space, as measured by an asymptotic observer 
living in an asymptotically flat space-time. We find that for Quantum Gravity 
(QG) it is energetically favourable to perform its quantum fluctuations not upon 
flat space-time but around a "gas" of wormholes of mass m p , the Planck mass 
(m p ~ 10 19 GeV) and average distance l p , the Planck length a p (a p ~ 10~ 33 cm). 
As a result, assuming such configuration to be a good approximation to the true 
Ground State of Quantum Gravity, space-time, the arena of physical reality, turns 
out to be well described by Wheeler's quantum foam and adequately modeled by 
a space-time lattice with lattice constant l p , the Planck lattice. 

PACS 04.60, ...(other pacs) 



1 Introduction 



Among the fundamental interactions of Nature, since the monumental contribu- 
tion of Albert Einstein, Gravity plays the central role of determining the structure 
of space-time, the arena of physical reality. As well known, in classical physics 
a world without matter, the Vacuum, has the simplest of all structures, it is flat 
(pseudoeuclidean) ; but in quantum physics? This is the focal question that has 
occupied the best theoretical minds since it became apparent, at the beginning 
of the 30's, that Quantum Field Theory (QFT') is the indispensable intellec- 
tual tool for discovering the extremely subtle ways in which the quantum world 
actually works. Thus the problem to solve was to find in some way or other 
the Ground State (GS) of Quantum Gravity (QG), which determines the dy- 
namical behaviour of any physical system, through the non-trivial structure that 
space-time acquires as a result of the quantum fluctuations that in such state 
the gravitational field, like all quantum fields, must experience. Of course this 
problem, at least in the non-perturbative regime, is a formidable one, and many 
physicists, J. A. Wheeler foremost among them, could but speculate about the 
ways in which the expected violent quantum fluctuations at the Planck distance 
l p (l p = 10~ 33 cm) could change the space-time structure of the Vacuum, from its 
classical, trivial (pseudoeuclidean) one. And Wheeler's conjecture, most imagi- 
native and intriguing, of a space-time foam vividly expressed the intuition that 
at the Planck distance the fluctuations of the true QG ground state would end 
up in submitting the classical continuum of events to a metamorphosis into an 
essentially discontinuous, discrete structure^ 

It is the purpose of this paper to give a detailed account of the results of an 
investigation on a possible QG ground state, which has been summarily reported 
in a recent letter [|IJ. The starting point of our attack is the realization that QG 
can be looked at as a non-abelian gauge theory whose gauge group is the Poincare 
group. Following the analysis performed by one of us (GP)|2j of another non- 
abelian gauge theory QCD (whose gauge group is SU C (3)), we decided to explore 
the possibility that the energy density (to be appropriately defined, see below) of 
the quantum fluctuations of the gravitational field around a non-trivial classical 
solution of the Einstein's field equations for the matterless world, could be lower 
than the energy of the perturbative ground state (PGS), which comprises the 
zero point fluctuations of the gravitational field's modes around flat space-time. 
Indeed in QCD it was found that the unstable modes (imaginary frequencies) 
of the gauge fields around the classical constant chromomagnetic field solution 
of the empty space Yang-Mills equations, in the average screen completely the 
classical chromomagnetic field, allowing the interaction energy between such field 



1 We should like to recall here that, based on Wheeler's idea, a successful research program 
was initiated a few years ago to explore the consequences of the Standard Model (SU C (3) <g> 
SUl(2) <g> Uy(1)) in a discrete space-time, conveniently modeled by a lattice of constant l p , the 
Planck lattice (PL). 
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and the short wave-length fluctuations of the quantized gauge field, to lower the 
energy density of such configuration below the PGS energy density. Thus we 
decided to try for QG the strategy that was successful in QCD, i.e. 

1. select a class of empty space classical solutions of Einstein's equations 
that is simple and manageable; 

2. evaluate the spectrum of the small amplitude fluctuations of the grav- 
itational field around such solutions; 

3. set up a variational calculation of the appropriately defined energy 
density in the selected background fields; 

4. study the possible screening by the unstable modes (if any) of the 
classical background fields. 

As for point (1) we have chosen the Schwarzschild's wormhole-solutions|3|], 
the simplest class of solutions of Einstein's equation after flat space-time. In 
order to achieve (2) the Regge- Wheeler ||] expansion has been systematically em- 
ployed, yielding a well defined set of unstable modes (for S-wave). This important 
result, already indicated in previous independent work ||, renders the develop- 
ment of the points (3) and (4) both relevant and meaningful, the former point 
yielding a lowering of the energy density due to the interaction of the short-wave 
length modes with the background gravitational field, the latter exhibiting the 
(approximate) cancellation of the independent components of the tensor of the 
Schwarzschild's wormholes by the S-wave unstable modes. As a result flat space- 
time, like the QCD perturbative ground state, becomes "essentially unstable", 
in the sense that upon it no stable quantum dynamics can be realized. On the 
other hand a well defined "gas" of wormholes appears as a very good candidate 
for the classical configuration around which the quantized modes of the gravita- 
tional field can stably fluctuate. But a discussion of the physical implications of 
our findings must await a more detailed description of our work, which we are 
now going to provide. 

2 The Schrodinger functional approach 

In order to develop a functional strategy aimed at determining the Ground State 
of Quantum Gravity, which parallels the approach developed for QCD [0, we 
must first identify an appropriate energy functional. In General Relativity this 
is a non-trivial problem for, as is well known, in the canonical quantization pro- 
cedure, first envisaged by Dirac [H] and Arnowitt, Deser and Misner (ADM) pj, 
due to general covariance the local Hamiltonian is constrained to annihilate the 
physical ground state, a fact that in the Schrodinger functional approach is ex- 
pressed by the celebrated Wheeler-DeWitt equation ||. However we note that 
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the problem we wish to solve concerns the minimization of the total energy of 
an "open space", in which there exists a background metric field that becomes 
"asymptotically flat", i.e. that for spatial infinity (r = |af| — ► oo) behaves as 



9^^vSP + 0(-) t (1) 



r 



where rfffi is the Minkowski metric. In the conventional canonical formulation, 
space-time is foliated into spacelike slices E with constant values of the time 
parameter t; the asymptotic condition (j^) determines the asymptotic behaviour 
of the canonical variables: the spatial 3- metric gij on E, the conjugate momenta 
ir ij , the "lapse-function" N and the "shift- vector" iV; |7J] as: 

9ij — > Si.+Oir- 1 ) 

n ij — ► 0(r- 2 ) 

N — ► 1 

N t — > 0. (2) 

Let us consider the ADM-energy [[I]], which in cartesian coordinates is given by 
(<9E is the boundary of E , ", k" denotes partial derivative with respect to x k , 
and G = l 2 is the Newton constant) 

Eadm = / dS k S ij (g iktj - g ijtk ). (3) 

167TG JdY, 

We should like to point out that Eadm is just the energy that an asymptotic ob- 
server attributes to a space region E whose time foliation he is keeping anchored 
to his (asymptotically) flat metric: the boundary conditions select a privileged 
reference frame (up to Lorentz transformations) that implicitly defines the phys- 
ically relevant energy. In particular, the asymptotic condition on iV fixes the 
"boundary time" unequivocally: the asymptotic observer is the only possessor of 
an idealized clock that allows him to describe quantities associated to the whole 
physical system without introducing material clocks (i.e. auxiliary fields). In 
this sense, he is also the only one that can really be termed as an idealized, 
non-interfering "observer" capable of describing geometry at the quantum level 
in terms of evolution, not merely in terms of correlation between variables, thus 
giving a full meaning to the expression "quantum geometrodynamics" 0. 

At the classical level, the definition of Eadm fixes Minkowski geometry as the 



zero point of the total energy; the proof of its positivity |T0[ can be looked at as 
the statement that flat space-time is the (unique) vacuum of General Relativity. 
This explains why the first steps towards the quantization of the theory have been 
based on a perturbative approach on the flat background, with the fluctuating 
self-interacting field interpreted in the conventional particle view as creating and 
annihilating gravitons, which propagate in pseudo-euclidean space: in this sense, 
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we call flat space-time filled with gravitons performing zero-point fluctuations the 
"perturbative ground state" (PGS). 

We now know that the perturbative approach was doomed to fail: the non- 
renormalizability of the theory does not allow to make any meaningful and pre- 
dictive perturbative expansion. On the other hand, analyzing the theory to the 
lowest non-trivial order around a curved background may give us important indi- 
cations of how the deadly "impasse" of the perturbative approach may be finally 
overcome and give back to the simplest form of QG its status and role of a "bona 
fide" Quantum Field Theory. 

We thus study the quantum fluctuations of the gravitational field upon a 
generic, asymptotically flat stationary background geometry, solution of the source- 
less Einstein's equations; in particular, for the background metric we can choose 
a foliation orthogonal to Killing timelike vectors and put it in static form. On a 
given slice E, the 3-metric is thus given by 

gij{x) =rnj(x)+hij(x), (4) 

where rjij(x) is the spatial background metric and hij(x) the fluctuation to be 
quantized (i G E). We can now expand the total ( the sum of the Hamiltonian 
and Eadm) energy E of space in powers of the fluctuations hij(x) (the number 
(n) denotes the order of the expansion): 

E = Ef DM + E% M + £ / d 3 x(NH^ + NtfM). (5) 

n>l JS 

7i and 7i l are the super-hamiltonian and super- momentum, as defined by ADM 
0: 



T = 167rGg mi 7r ij TT kl 



n = t + v 

"J 

where g is the determinant of 3-metric and R the corresponding curvature scalar, 
with the "supermetric" Giju given by 

Qiju = -g~ 1/2 {gikgji + gugjk - gijgn), (7) 

while (; is the covariant derivative with respect to #y) 

TC = -27r^.. (8) 

We note that, a priori, also iV and Ni can be expanded in (H); while the back- 
ground terms iV^ ) and are fixed functions, subject to the asymptotic condi- 
tions (El), the higher order terms represent true fluctuations in the lapse function 
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and the shift vector: variations of N and N{ yield, at the classical level, the 
constraints: 



n = n(n) = 

n>l 

H i = Y / n i {n) = o, (9) 



n>l 



leaving in the expression of the classical energy only the E^ DM term. 

Since 77^, A r< -°- ) and form together a solution rj^ of matterless Einstein's 
equations, the linear term in the canonical Lagrangian density must be a total 
divergence. In our case, where rj^ u is static, this is a purely spatial divergence 
and, keeping the asymptotic flatness of background in mind, it must necessarily 
coincide with E^ DM . Thus, we have: 

e% m = ^ / gs ds k (h kjij - h jj;k ) = - jf d 3 x(N^n^ + nPw^i (10) 

that, together with the constraint (§) allows us to rewrite the energy (H) as 

E = Ef DM + J2 I d 3 x(N^H {n) + iVf } ?f (">). (11) 

n>2 ^ 

Note the survival of only the classical (O)-order terms in iV and N{. 

The quantization of the theory promotes the canonical pair hij(x),ii kl (y) on 
X to operators obeying the commutation rules: 

[hv(x),iP(y)]=ih6l k $6®(x-y), (12) 

and acting on a Hilbert space of functionals \I/ that are annihilated by the 

constraints (Q). The evolution of the physical states ^ is governed by the 
"Schrodinger equation" 

ihd t ^> = E [hij,ir kl ] (13) 

where the Hamiltonian operator is given by (|TT|). Note that this is just the 
description of the quantum dynamics made by the asymptotic observer at infinity. 

We point out that our definition of the Hilbert space is truly consistent within 
our restriction of phase space to two pairs of canonical operators, obtained by 
the gauge conditions (that respect ([5])). This does not mean a loss of invariance 
(and of physical reality) at all: despite its look, the total energy (|TT|) is nothing 
but a rearrangement of the invariant ADM-energy. 

We should also be aware that the configuration space representation of the 
canonical operators 

hij(x) — ► h i:j (x) x 

*" ix) - - !fi «4r (14 > 
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acting on state functionals 

tt — tf[/Hu](i) (15) 

is not easily manageable beyond the 1-loop level, where connected ghost terms 
appear. Beside that, for n > 3 the expansion of the Hamiltonian operator (O) 
contains products of conjugate operators, thus posing an ordering problem. These 
problems are related to the bad ultraviolet divergences that would still yield a 
non-renormalizable behaviour, the possible solution of which emerges from the 
results of Section 6, which show that the structure of the vacuum is, with good 
probability, essentially discontinuous at the Planck scale l p . Thus, in the rest 
of our analysis the QFT we shall work with will be cut-off at the Planck scale, 
having clearly in mind that our results will only be meaningful if consistent with 
this fundamental assumption (see Section 6). 

As for the constraints (|Sp, the problems are easier to solve. In fact we first 
notice that at the lowest order, the Hamiltonian operator retains only quadratic 
terms in the fields, on which we have to impose consistently first order constraints, 
that do not annihilate the quantum energy. The following terms in the expansions 
(HD can be enforced through a systematic correction of the state functional \1/ that 
readapts non-physical degrees of freedom order by order, thus not affecting the 
dynamics based on the degrees of freedom (two for each space point) isolated at 
the lowest level. @ Thus in spite of the problems typical of Quantum Gravity, the 
parallelism with the situation in QCD [[| is fully regained. 

According to our fundamental assumption to cut the theory at the Planck 
scale, we shall perform a 1-loop calculation, with the Hamiltonian operator trun- 
cated at n = 2. We simply adopt the representation ( |bj|) and (p3[), thus obtaining 
the Schrodinger equation 

ihd t ^[h t ,j\(t) = H[h tJ) -ih—-Mh ij ](t), (16) 

Ollkl 

where 

H = M+ [ SxN^n^HhiM^-th „ 6 , A (17) 

JS 0llij{X) 

with ty[hij](t) annihilated by the first order constraints 

n^\h tJ (x),-zh JJ ^ ) Mh lJ }(t) = o 

W (1) (M*), -ih-jj^Whijlit) = 0. (18) 



Setting, as usual, 



Vfc]® = e-™'**^), (19) 



2 See for example the procedure followed in Appendix B of ref 
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the problem can be reduced to the eigenvalue equation 

H[h ih -ih-^-Mhij] = EV[hij]. (20) 

We can now investigate the ground state of the theory. Instead of solving directly 
the eigenvalue equation, we look for the minimization of the expectation value of 
H on a class of gaussian wave-functionals: 

E m = J [vh]^*[hij\H 

If the background is stable under the action of quantum fluctuations, at the 
1-loop level this result coincides with the solution of fl20|); if, on the contrary, 
simple minimization leads to an imaginary part in E^ , then we have discovered 
an unstable configuration, whose physical meaning must be investigated. We 
demonstrate in the next section that the latter case occurs when fluctuate 
around the "wormhole solution" 77^ discovered by Schwarzschild in 1916 0, whose 
line elements in polar coordinates are given by (2GM < r < +00) 

ds 2 = - T ~ 2GM dt 2 + dr 2 + r\dd 2 + sin 2 6d<p 2 ) (22) 

r r — 2GM 

and depend on the single parameter M, the ADM-mass, such that 

Ef DM = M. (23) 



(21) 



3 Quantum Fluctuations on a Schwarzschild 
Background 

We shall now address the problem to evaluate the expectation value on a gaussian 



trial functional of the Hamiltonian ( |I7D where, according to our fundamental 



hypothesis (to be checked for consistency at the end of the calculation), we keep 
only the quadratic terms in the field quantum fluctuations hij. This truncation 
corresponds to the one-loop approximation. From a classical standpoint this 
amounts to a calculation of the energy carried by the quantized gravitational 
waves propagating on a fixed background, in the weak field approximation. 

In our analysis we shall follow closely the steps of the ref . , where a similar 
calculation was carried out for a Yang-Mills theory. We begin by constructing 
the Hilbert space of the states of the gravitational field, introducing the following 
scalar product: 

((*|*')> = / (Vh)A FP ^[h tJ ]^'[h tJ }, (24) 

where (T>h) denotes the measure of the functional space and App represents 
the Fadeev-Popov determinant, depending on the gauge adopted, necessary to 
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recuperate the gauge-invariance, i.e. the general covariance of QG. The Hilbert 
space will thus be the space of the state- vectors *S?[hij](t) that with the metric ([24]) 
are normalizable. We note that for an infinitesimal coordinate transformation the 
quantum field hij gets transformed as: 

h i:j (x) — > h i:j (x) - &\j(x) - £j\i(x), (25) 

just like a weak classical field. And in our approximation, being the gauge- 
conditions (H) linear in the field h^, the determinant A FP does not depend on 
and can be therefore neglected. For a generic operator 0[hij, n ], the expectation 
value on a state ty[hij)(t) can be defined as: 

((0)) * (t) = ~W\W ( } 

Let us consider now a hypersurface £ at a fixed time t. We wish to com- 
pute the expectation value of the (truncated) Hamiltonian on the gaussian trial 
functional: 

* G [^-] = exp -\f JHjWrW&flhuOf) (27) 

4 Jx,y 

where J x = / s <i 3 x?7 1//2 (x) and r/^ (x) = r ^ 2m (m = MG is the one half Schwarzschild 

radius). In order to get a normalizable tyc[hij] we require that F v (x, y) be real 
and positive, symmetric under the exchanges i <-» j , k <-> / , ij,x <r-> kl,y. 
The second order Hamiltonian density is given by 

ft( 2 )(£) = T (2) (x) + 2) (£) (28) 

with 

= -167rGtT 1/2 faV - U^ 1 ) _f (29) 

2 dhij^xjdhkilx) 

V {2) = -^v 1/2 (R (2) + \h k k R {1) ). (30) 

lb7rG 2 

We observe that for a Gaussian wave-functional the expectation value of two 
fields is given by: 

({hij{x)h kl {y)))<s> = G ijk i(x,y) (31) 
where Gijki(x,y) satisfies the relationship: 

jf r^ mn (x, z)G mnkl (z, y) = r^WlA^^ - V)- ( 32 ) 

In this way one gets for the expectation value of the Hamiltonianll 

a ■ 1 1 

-(r^.(f,f) - -rv. (£,£)) + * 

4 2 a 



((H))* = M + h j N{x)[^ ^x) - -rVi&a)) + -0^(x)G^-(x,x)], 



(33) 



3 In order to clearly separate the classical from the quantum (one-loop) contributions for the 
rest of this Section we shall keep the Planck constant h, instead of putting it equal to one, as 
done in the natural unit system. 
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where 

a = IQirGk = lQnl 2 P , (34) 
and the differential operator 0^ kl (x) is defined by 

NV® = - [ N(x)h l3 (x)& k \x)h kl (x), (35) 



a Jx 

and represents the "potential" contribution to the quadratic Hamiltonian (|17|). 
In order to guarantee the general covariance of our calculation, it is necessary to 
impose on the physical state the quantum constraints, which in our approximation 
are 

H^^ = - h ij {ij - /// nl/ = (36) 

7^(1)^ = ^ = o (37) 

which are obeyed provided, 

V,(f)P^(f,f) = (38) 

where Vj(af) denotes the covariant derivative with respect to the background 
field; and, fixing the lapse function and the shift vector as 



N = = W 1 - ^ , N l = N*® = 0, (39) 

we have 

V,(|) = 0, h k k = 0. (40) 

By consistency with the trace condition h\ = we must also impose 7i k ^/ = 0, 
which yields the further constraint: 

T t M (x 1 y) = T^ k k (x,y) = 0. (41) 

The elements of the Hilbert space of the physical modes of the gravitational 
field hij are thus the symmetric tensors of rank 2 0^, defined in S, obeying the 
boundary conditions @ and the gauge conditions: 

Vi(^) =0,H = 0, (42) 
normalizable with respect to the scalar product: 

($!$') = / iV- 1 ^* (43) 

Jx 
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In this way we may construct in our Hilbert space a complete orthogonal system, 
by making use of the spectral decomposition of the operator O y (x) or, better,, 
of the operator Q^ kl defined as: 

NV m = I [ N(x)h l3 (x)d^ M (x)h kl (x) = I N-\x)h l3 {x)Q i ^\x)h kl {x) ) 
a Jx Jx 

(44) 

where in order to go from O to Q a total divergence has been added to the 
integrand, without changing the "potential" contribution to the Hamiltonian. 
Thus the operator Q becomes in our Hilbert space a self-adjoint hermitian 
operator, whose eigenf unctions (p denotes a complete set of indices), 

= A(p)$i p) , (45) 

build the sought complete orthonormal basis. In this basis the "propagator" 
Gijki(x,y) has the simple form; 

G m (x,y) = £ 2^y*lf GO, (46) 

where f(p) denotes a set of variational parameters to be determined by the min- 
imization of the expectation value fl33|). 
From (|32|) and (^) we obtain, 

T ijkl (x, y) = £ 2f(p)N-\x )N-\m\;i{x)<Pl ) {y). (47) 



We may now easily compute the expectation value (p3|) , and obtain: 

((H)h = M +\ E(«/(p) + T^y)- ( 48 ) 

and minimizing with respect to the variational function f(p), i.e. imposing 

* ((H)h = 0, (49) 



sf(p) 



we readily get 



f(p) = a~ 1 ^X{p). (50) 
which inserted in (|48|) finally yields: 



E = M + n^y/\ip). (51) 
p 

All the above makes sense if and only if 

Hp) > o, (52) 
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i.e. the eigenvalue of the "potential" operator Q are positive definite. If, instead, 
for some p X(p) < the one-loop approximation, yielding imaginary contributions 
to the energy of the ground state, breaks down, showing that the Perturbative 
Ground State (PGS, M — > 0) is essentially unstable. 

This is precisely the situation found in the study of SU(n) Yang-Mills theories^] 
where, going beyond the one-loop approximation, one could easily check that the 
modes belonging to the sector where A(p) < did not contribute to the energy 
of the state terms of 0(/i)i, but rather of 0(1), just like the classical term M. 
This "promotion" of a quantum 0(h) contribution to a classical 0(1) one, can 
be understood when we realize that the amplitude f(p)~* of the modes with 
A(p) < 0(the "unstable modes") [see Eqs. (|3TD and (|||)] is only prevented from 
becoming infinite by the neglected positive terms of 0(h 2 ). In this way ctf(p) 
becomes O(^) and the negative contribution from the "unstable modes" is just 
classical, i.e. 0(1). 

In the calculation of Ref. @ one could explicitly prove that this "promoted" 
quantum contribution completely screens the classical positive term (such as, in 
our case M), thus realizing a "vacuum" state whose energy density is way that 
of below the PGS, which as a result becomes unstable at all space-time scales. 
In the case of QG the problem of going beyond the one-loop approximation is 
formidable, utterly beyond our present means of analysis, however, as shall be 
discussed below, to figure out the contributions to the energy of the trial states 
of possible "unstable modes" with X(p) < appears reasonably doable. 

In order to precede further we must first compute the operator Q l ^ kl (x) and 
then diagonalize it, which we shall do next. Defining = — rfi, we have: 

W = -h ij + h\h lj - h}h l m h mj + ... (53) 

and for the perturbative calculations it is useful to introduce: 

oa _ -pa _ -pa (s) (xa\ 

^mn mn mn \ ) 

which turns out to be tensor. 

Observing that (";" denotes the full covariant derivative, while "|" is the 
covariant derivative with respect to the Schwarzschild background) 

hja-,i = hj k \i — S l im hij — S\jhik, (55) 

we may write the Riemann and the Ricci tensors as: 



pa p a ( s ) I ca ca i qa Cm Ca cm 

^bcd — n-bcd °bd\c °bc\d ' °mc D bd D dm°bc 

p p ( s ) _i_ ca ca j_ ca cm ca al 

t^bc — J^bc ~r °bc\a °ab\c D am°bc D lb°ac 



(56) 



4 Like it happens, according to (|5l|), to the modes belonging to the "stable sector", for which 
A(p) > 0. 
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furthermore: 



Sjk - ^9 %k {hjr\k + Kk\j - h jk \r)- 



(57) 



We may now compute the curvature scalar R = g bc R hc , which for convenience we 
decompose as the sum of three terms: 

A EE g b <Rl s J; B ee - g^ C ee g bc (S l lm S™ - S l mb S™) . 

In this way to first order in hij we have: 

A W = ^(-h + 3h{) (58) 

while to second order we obtain: 

A (2) = ™ 3 (h«K-3h}h\). (59) 

As for B we have 



d(1) _ J.mqa(X) _ J-mq^iX) 



iab _ t a 

It \hr, lt\ „ 



(60) 



and 



B {2) = { ~ [h ar Kj s ] + \h^ a h\ r - h lm h« mlal } + {h lm h lm - 2hA^} 



-3A& + h*»hj? a + \h\h x £ 



respectively. While for the third term one has: 

C {2) = \h\ b h b l - l -h\ b h} b - \{2h al \ b h aW - h am \ b h am \ b ). 



(61) 



(62) 



By summing the different terms we obtain the following expansion of the scalar 
curvature: 

R {0) = (63) 



= {h ab lba } 



m u U a 



(64) 



R^ 



h ar hj s 



\a 



(65) 
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But we are not done yet, we must expand the square root of the determinant of 
the metric J \\g\\- Due to = we may stop at the first order, obtaining: 



\w\ = 1+ ^ hl+0{h2 ' ) 

In this way for we get the following expansion: 



(66) 



V {2) 



IQttG 



2r 3 



+3^h\h - (h ar hj s ) ]a + h ar ]a h lr - h lm h a m]al + h lm h llm 



+h lm h 



r 

Imu \ a i h lm h ' r L h al h b 
'Ira a "T a' 1 \r'Hm o \b a a\l 



(67) 



We are now ready to compute the operator Q^ kl . Let us fix the gauge as: h — 
and (Jjf*) — 0, the latter condition by use of the former giving 



We also note that given a vector T m , one has 



N[T m ] {n 



- / T 1 



1 



rn 



Thus we obtain for the integral of the second order potential 



(68) 



(69) 



| 6 m 



2 a l|«^2 



c D 6 



r 3 (r — 2m) 
Recalling Eq.(pH]), Eq.(|45|) becomes: 



(70) 



iV 5 



r 3 (r — 2m) 



A 1 h b) 



(71) 



where we have added the terms in v 1 e r in order to remain inside the Hilbert 
space of the tensors 4> l \ which obey the conditions [see Eq.fl42|)1 



li 



'Hi- 



(72) 
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By taking the trace of eq.([71]) with 77^ we get immediately that 



T 



-v b n 2 (rW 10m2 ^ ^ 

3 16 1 11 3r 3 r-2m 



(73) 



By consistency, we must require that the terms in v j do not contribute to the 
potential energy, this obviously implies that 



\j 



0. 



(74) 



which requires that for r — > 00 the following conditions be satisfied: 

vi -> 0(r- x ) 
w 2 ,w 3 -> 0(1), 

and 

lim r 2 / dQ<f)\vi = lim r 2 / dQ<j)\vi = 0. 



(75) 
(76) 



4 The eigenvalues and eigenmodes of the sec- 
ond order potential 

Due to the spherical symmetry of the Schwarzschild background a particularly 
suitable method to obtain the solutions of the eigenvalue problem posed by 
Eq. (|7T|) is the one devised by T. Regge and J. A. Wheeler j|, for the study of the 
small (classical ) fluctuations around the Schwarzschild solution. By making use 
of this method the eigenfunction are separated in two classes (see appendix B), 
the "even" solutions with parity (—1)', equal to the parity of spherical harmonics 
^m(#,0), and the "odd" solutions with opposite parity.For the "even" solutions, 
if we set 

' <fi = H{r)Y lm (9,<p) 
<Pl={G 1 {r) + G 2 {r)df)Y lm {e^) 
<%= {Gi(r) + ^}dl + G 2 (r) cot Ode) Y lm (0^) 
<P\ = K{r)d e Y lm {6^) 

. 4>l=^{de- cot 9)8^9, if), 

the eigenvalue equations will turn out to be, as we shall see in a moment, com- 
pletely factorized. As for the vector Vi ([H]), factorization is achieved if we set, 



v x = u{r)Y lm {9,ip) 

v 2 = V(r)d e Y lTn (9,<p) 
v z = V{r)d 9 Y lrn {9,ip) 



(7f 
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Substituting ([77]) and ( [78] ) in ( |7I| ) (see appendix C) we obtain a system of differ- 
ential equations for the radial functions only: 



' N 2 { 



r—2m Q2 _|_ 2 r ~2m 



a 



2 6r 2 -27mr+23m 2 



#(r) 



+ 4^Z(Z + l)K{r) + 2^(2G?i(r) - /(/ + l)G 2 (r))} 



-- -H(r)A\ + I 



r—2m 



d r U(r) 



r—3m 



U(r) + 



V(r) 1(1+1) 



N 2 { 

= -A\G 2 (r) + 2 -^ + N 2l -^G 2 (r) 

iV 2 { [^dl + 2^9 r - d(r) + ;|G 2 (r)/(Z + 1) 



(79) 



+ 



2 3r 2 — 12mr+7m 2 

3 r 3 (r-2m) 



ff(r)} 



= -4AG? 1 (r)-| 

+jV 2lM Gi(r) 



r— 2m 



d r U(r) 



r—3m 



U(r)-mi(l + l) 



N 2 { 



r—2m f\2 _i_ 4i — 7m Q 2r 2 — mr-10+lOrre 2 

r ^ T r 2 r 3(r-2m) 



AT(r) 



+ 



2r— 3m 



^-^yF(r) - £ [d(r) + (1 - 1(1 + 1)) G a (r) 
-A\K(r) + i ([/(r) + d r V(r) - f V(r)) + iV 2 ^±Ii^(r). 



And by substituting the same expansions in the constraint Eq.(^), we obtain 

2G x (r)-G 2 (r)l(l + l)+H(r) = 
3H(r) m 



d r H{r) + 



-#(r) - AT(r)Z(Z + 1) =0 



r(r — 2m) 

d r K(r)r(r — 2m) + i^(r)(4r — 8m) 
+G 1 (r) + G 2 (r) [1 - 1(1 + 1)] = 0. 

The boundary conditions (|75]) for r — > oo become: 

' U(r)^0(h 



0) 
1) 



(82) 



J3) 



I V(r) -0(1) 

while ( f76| ) imply: 

{r 2 [if (r) • C/(r) + 1(1 + l)K(r) V»] } 

= {r 2 [#(r) • £/(r) + 1(1 + l)A"(r) V(r)] } 



I r=2m 



0. 



»4) 
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Proceeding in a completely analogous way, for the "odd" solutions one gets: 
'01 = 

02 = F 2 (r)^[d -cote]d v Y lm (e,ip) 



01 = 



-01 



(85) 



0? = SS^w^) 



j.3 _ 1 J2(r) 
l ^2 — 2 sine 



^aj + cotfl-^liu^) 



and 



V! = 

v 2 = D(r) (-^dyY^v) 
v 3 = D(r) smed Y lm (e,<p). 



(86) 



And substituting these last expressions in the eigenvalue equations, we get 
the following relationships: 



N 2 



r-2m 2 6r-llm a 1(1 + 1) , 2 



r 

r — m 



d r + v - y + — (-3(r - 2m) 



,y*2 .y*3 



r — 2m 



-m 



F 1 (r) = 4AF 1 (r) + ^(d r --)£>(r) 



(87) 



and 



N 



r-2m a2 r-3m a Z(Z + 1) 2\ 



+2- 



r 

2r — 3m 



^•?2 ^<2 



i^i(r) 



4AF 2 (r) - iiXr) 



Noting that the fluctuations' tensor automatically satisfies the traceless con- 
dition, we find that the gauge conditions yield 



(r - 2m)(r<9 r + A)F 1 (r) - ( 1 - ^tli ) p 2 (r) = 



(89) 



as the only constraint. 
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The conditions ([T5|) reduce to D{r) — > 0(1), while ( |T6[ ) yield 



/(/ + l)r 2 F!(r)D(r) 



|r=2m 



/(/ + l)r 2 F 1 (r)D(r) 



0. 



(90) 



An important observation is that for S-waves (I = 0) the gauge conditions are 
sufficient to determine the form of the solution. In the other cases we shall solve 
the eigenvalue problem in the WKB approximation. 

Setting now I — in our equation we realize at once that the "odd" part 
vanishes identically, while the first of the (|7^) has a particularly simple form, de- 
pending on H{r) and G{r) only, the latter obeying through (^0|) H{r) = —2G\{r). 
The fluctuations' tensor acquires thus the following simple form: 



cf ) l^6^) = -^{z5 a l 5l-5l) l -H 
where H(r) is determined via eq. (jSlp 
d r H(r) + H(r 

whose solution is 



3 1 /l 
r + 2 



(91) 



r 



1 

r — 2m 



0. 



(92) 



where the constant A is determined by the normalization condition. 

A 2 



1 



J d 3 x\\rf 



NAk2 



16m 3 



(93) 



which yields A = 4m 3//2 , so that 

H(r) = 4m 3/2 i 



2m 1 



3 ' 



(94) 



In Appendix D, we develop the calculation for the eigenfunctions' normalization 
in general. 

It remains to be verified whether there exists an eigenvalue corresponding to 
(|94"D, in other words, we are looking for a value A for which (0) is a solution of 
the eigenvalue equations satisfying the boundary conditions ( |75| ) and ([76]) . On 
substituting (|94|) in the first of (|79|), and observing that 



d r H(r 



3r — 7m 



H(r) 



d 2 H(r) 



r(r — 2m) 

12r 2 — 56mr + 63m 2 
r 2 (r — 2m) 2 



H(r) 
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and 



we find: 



G 1 (r) = --H{r) 



4 / r — 2m d r — 3m \ 
~ 1 — u ( r ) 



r dr r 2 



whose solution is: 

U(r) 



f _ m(12r-35m)\ 3/2 j_ 



I ^,2 2r^ 4 I 



(95) 



r — 2m 



C being an integration constant. Due to 
fixes the eigenvalue A. Indeed from 



51) C must vanish, while 



(96) 



finally 



lim r H(r)U(r) 

r— >2m 



one gets for the eigenvalue A: 



A 



1 



64m 2 64(GM) 5 



(97) 



This finding of a negative eigenvalue A, i.e. of an "unstable mode" around a 
"wormhole" , in view of the discussion in section 3, is a most significant result of 
our analysis that strengthens the analogy of the present calculation with the QCD 
one0. We must, however, point out that (P7|) is not at all unexpected, since in 
ref.0 a similar finding has been reported in a somewhat different context. 

As for the remaining stable modes we shall solve the eigenvalue problem in an 
approximate way, by the WKB method. We shall thus obtain a good description 
in the semiclassical region. Let us then write our solution in the form 



( H{r) 
K(r) 
Gi(r) 
G a (r) 
U(r) 
\ V(r) ) 



( h(r) ^ 

k(r) 

9i(r) 

9*{r) 
u(r) 

\ v ( r ) J 
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1 ho(r) 


+ (-) hAr) 


+ 




Mr) 


+ (!) Ur) 


+ 


(j)Mr) + --- 


9i ( r ) 


+ {j)9n(r] 


+ 


(?)V a (r) + --- 




+ (f)^(r; 


+ 


(!)% 22 (r) + --- 


u (r) 


+ (!) 


+ 


(f) 2 M2 (r) + --- 




+ (f) Mr) 


+ 


(!) 2 - 2 w + --- J 



• e 



i(so(r)+f S i(r)+(f) 2 S2 (rH 



(9? 



In order to simplify the analysis we enclose our system in a spherical box of 
radius R^> m, and impose the condition that the field vanishes at the boundary, 
in this way the energy spectrum becomes discrete. By using the gauge constraints 

h 

and equating the terms of the same order — we have 



( H{r) 
K(r) 
Gx{r) 
C 2 (r) 
U(r) 

{ v(r) ) 



( 



(!) 



i(i+i) 

2 



s (r) N 2 



iff Mr) 

£2flM ^ + (!) 2 ^ 2 ( 



^9Mr) 



h 2 (r) 



g2 (r) + (!) g 2l (r) + (f ) <? 22 (r) 
«o(r) + (!) Ui(r) + (f) 2 M2 (r) 
^o(r) + (f)%i(r) + (f )% 2 (r) 



e l(^(0+(!) S iW+(?) 2 S 2W 

with the conditions (/ = 4-f(r)): 



(99) 



fc 2 (r) 



s (r) 
1 - 



so(r) 

Hi(r) + ^i(r)+ 52l (r) (l 

32 



+ Mi(r) 



(100) 



i(i+i) 

2 



1 



so(r) N 2 r 2 ' 

Inserting these expressions in the eigenvalue equations (|79f) , we can check that 
the following is a solution: 

/ 



/ H{r) ^ 
AT(r) 
Gi(r) 
G 2 (r) 
U(r) 
{ V(r) ) 



±(!) '(i + i) 



i 



± 2 



/(*+!) 
2 



2v / Ar 



2 



t(!) + 
i 






(f) : 

1 - 



JV 2 
4Ar 2 

_ 1(1+1) 
2 



\ 



1 

4Ar i 



1(1+1) 
2 



N 2 
4Ar 2 
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2^ +pVXr*+a(l)+0(n) 



N--e 
r 



(101) 



where c is an integration constant, depending in general on I, and r* = r + 
2mlog^2^; — l)- The constant a(l) is crucial for determining the eigenvalues. As 
illustrated in Appendix E, a{l) can be determined by comparing ( |101| ) with the 
exact solutions for the flat case M = 0, obtaining 



a(l) = + 



71 



We shall now determine the regime of validity of our approximation. We observe 
that ( |101[ ) has been obtained by expanding the eigenvalue equations up to the 
order (Wj . By considering the next order one gets: 



S2(r 



r 2m r 



2m 



(102) 



so that ft|s 2 (f)| <C 1 if Z^jf- ^ 1, implying that, as expected, the approximation 
is good in the very high energy region. Before imposing the boundary conditions, 
let us analyse the behaviour of the "odd" part. 
Looking for a solution of the form: 



( Fi(r) \ ( L{r) \ 

F 2 {r) = M{r) 

V D ( r ) ) V d ( r ) ) 

( L (r) + f Li(r) + (f )' L 2 (r) ^ 

M (r) + fM 1 (r) + (f) 2 M 2 (r) 
^ rf ( r ) + frf 1 (r) + (f) 2 d 2 (r) ) 

. e ±A (r)+A 1 (r)+f A 2 (r) 

and setting M = 1 and Mj = 0(j > 1), we easily find the solution: 

(F 1 (r)\ /±(f) 

F 2 (r) = 
V D(r) ) { 

. Ar 2C e± ±2v / Ar*-ia+l)|+0(ft)_ 

r 



(103) 



) 




i (n) 2 


[l '('+1)1 




2 


2^/Xr 2 V * J 


2 



AT 2 \ 
4Ar 3 \ 



(104) 



Let us first consider the flat case M = 0. By imposing the vanishing of the 
solutions, both "even" and "odd" for r = R, we obtain the spectrum: 



-R-l{l + l)\=n^ 
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or 



jKi = 7^(n + l + l)-. (105) 



2fT ' 2 

For M^O, one gets instead: 

which clearly shows that, due to that presence of the wormhole, the gravitons' 
spectrum is redshifted by the quantity jpr = - ^ % -v . In the next Section we 

shall analyse the consequences for the energy of our state of the results obtained 
so far. 



5 The energy of the quantized gravitational 
field around an ensemble of wormholes 

We have just seen that an external observer who looks into a finite spherical box 
(of radius R) , centered around a wormhole, perceives that the gravitons contained 
in it are redshifted with respect to the gravitons in the absence of the wormhole. 
We wish now to determine the over all energy shift 

A£7(M) = E(M) - £7(0) (107) 

between the gravitational quantum state in the box containing a wormhole and 
the one with no wormhole. In Section 2 we have seen that 



E(M) = M+ £ K,,i,m (108) 

n,l,m,tj 

where I, m are the angular quantum numbers, n is the radial one and a = =F for 
"even" or "odd" solutions. Please note that we have come back to the natural 
units where % = c = 1. In flat space-time the energy is 



E °= E ( 109 ) 

n,l,m, a 

Thus, from ( |105| ) and ( |106| ) we compute 

A£7(M) = E(M) - £7(0) ~ M - ^ m ^j~ E (n + l + 1). 

n,l,m,cr 

We must now sum over the quantum numbers, we get 

(n + l + l) = J2 2 ( 2/ + l)(n + l+ 1). (HO) 

n,Lm,cr nl 
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and introducing the momentum cut-off a! we have, 



From: 



AT 4 1 7 N 
E 2(2/ + l)(/ + n + l) = ^-^iV 3 + ^ 

n+l+l<N Z 10 J 

~ ^ = ^A*. (112) 



2 (ft 



7T 



we may finally write: 



64A 4 /? 2 /? 
^(M) = M--^GMlo g — , (113) 



or equivalently for the energy density: 

A£(M) M _36A\CM _R_ 

which clearly shows that the "stable" gravitons' modes on the Schwarzschild 
background, due to the gravitational red-shift, have a zero-point energy smaller 
than in flat space-time, leading to the negative term in flllol ). However we shall 
not analyse ( |113| ) any further, for this expression is deficient in two important 
aspects: (i) it does not incorporate the (negative) classical contribution from the 
"unstable" mode fl96|); (ii) it only describes a single, isolated "wormhole" which, 
in view of the large energy gains involved, if certainly does not correctly describe 
what we are after: a possible QG Ground State. In any event, the importance 
of (|113|) lies in its exposing unequivocally the quantum mechanical instability of 
the classical Ground State for the matterless universe: the (pseudo) euclidean 
space-time. 



6 The mult i- wormhole state 

In order to have a realistic model for the Ground State of QG, as we have just 
argued, we must seek a more realistic classical solution of matterless Einstein's 
gravity. It is clear from what we have learned so far that the natural candidate 
is a "condensed" system of WH's, either gaseous or crystalline, of ADM mass M 
and (average) density (^) 3 , a being the interwormhole average distance. Even 
though no explicit solution, that we know of, has been given for the general 

5 The problem of the cut-off, as already discussed, is a very important one. For consistency 
of our approach requires G(A)A 2 < 1, where G(A) is the value of the Newton constant at the 
cut-off A. 
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case, one such classical state must certainly exist, for by increasing a one is 
led to the well-known Schwarzschild solution, and special configurations of the n- 



wormhole system have been studied and characterized in ref. [11]. From this study 
one confirms the intuitive expectation that the ADM energy of "wormholes" of 
mass M, for separations a ^> 2GM (the Schwarzschild radius), coincides with 
the classical energy of a system of n matter points interacting via the two-body 
Newtonian potential U = Classically this configuration is unstable, since 

under the action of the Newtonian attraction the "wormholes" will not keep their 
(average) mutual distances a but will tend to collapse to a single "wormhole" 
configuration, encompassing the given space-time VE. Quantum -mechanically 
the situation is different, for the Heisenberg principle puts in general bounds 
upon the interwormhole distance a. Indeed, for a two-wormhole system, from the 
Schrodinger equation 



2 



u{d) = Eu{d), (115) 



where the reduced mass // = 4f, we obtain for the lowest lying state: 

GM 2 . . 

E = ~^- (U6 > 

with the "Bohr-radius" 

°° = G~Jp' (117) 
Thus, provided the "Bohr radius" a Q exceeds the "coalescence" distance d c = 
4GM (twice the Schwarzschild radius) the quantum mechanical system of n 
wormholes is stable. This means that one must have: 

a>a Q > 4GM, (118) 

the maximum allowed density occurring when 

a = a Q = AGM, (119) 

i.e. for 

M=Q 4 G-l (120) 

Extending (|113| ) to the n "wormholes" case, and setting R — |, we obtain for 
the energy density difference (a > 4GM — ) 

AE M 



3 



V a 



GM 16 2 4 / a 



4a 7T 3 \AGM 



(121) 



6 Evidcntly part of the energy in the collapse must be radiated away, as the single wormhole 
mass in the volume V is proportional to nsM and not to nM. 
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Thus, the energy density gain is maximum when a is minimum, i.e. when, ac- 
cording to (|117|) 

a = a ° = G^" ^ 

On the other hand, the quantum stability of the n-wormhole configuration teaches 

us that the minimum value of the Bohr radius a Q must equal 4GM, implying 

that the maximum energy gain is achieved when the wormholes have, according 

i i 

to (|1 20|) , an ADM-mass M = f|J 4 G~^ = (j^ 4 m p , and an average distance 



°"=^=G) iG " i =G) ii - (i23) 

We have just seen that, even without taking into account the (negative) con- 
tribution to the energy density of the unstable modes, a system of "wormholes" 
of ADM mass M = 2~^m p (m p = Gs ~ 10 19 GeV, the Planck mass) at the av- 

7 

erage distance f |J 4 l p (l p = G^ ~ lCT 33 cm, the Planck length) realize, according 
to fljgjp , a large energy (density) gain with respect to the Perturbative Ground 
State (PGS) of Quantum Gravity. 

We shall now try to answer the important question of the so far neglected 
contribution of the "unstable" mode around a single wormhole. Let us consider 
around a wormhole of mass M, of Schwarzschild radius 2m = 2GM, a spherical 
region of radius R = Am (— > A > 2). The spherical region shell: 2m < r < Am 
is the region (modulo the quantum fluctuations of the wormhole), which the 
asymptotic observer can probe, and in particular measure its (average) Riemann 
tensor. 

The existence around the wormhole of an " unstable" solution of the grav- 
itational field (which will obviously generalize to the case of many wormholes), 
found in section 3, guarantees us that such a mode will contribute a classical 
(negative) term to the ADM energy of the quantized gravitational field. What 
kind and size of contribution? As argued in Section 3, without involving ourselves 
in the intractable expansion of the QG Hamiltonian to 3rd and higher orders, we 
may figure out the effect of the "unstable" mode in the following way. Let us 
call the normalized unstable mode, which according to the developments in 
section 4, can be expressed in terms of the function H(r) given ( |94D and the con- 
straint function U(r), appearing in fl96D . Its contribution hijix) to the quantized 
gravitational field is thus 

hi = A^, (124) 

where A is a real constant whose value shall be determined in such a way as to 
minimize the Riemann tensor associated with the metric field 

9a >ih ■ h l- ( 125 ) 
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Thus we examine the metric: 



<,u = M - —) -2A 



,-3 



1 - 



Ira 



g 22 = r 2 + AJl-^r 1 



933 = | r' + A^l-^r-'j (sinOf 

which yields the non-vanishing components of the 3-dimensional Riemann tensor 
are: 



R 121 2 = {l23m 2 A 2 r 3 + 12m 2 A\l T -±^ - IW m A 2 r i I ////•'" 



+14 r 



7 —r + 2m 



mA -8mA 3 rJ r _±l^ + 2 7 A 2 r 5 



4r 4 / r _±J^ 8 Ar+lQAm X ' 



-r+2m 



(126) 



p 3 

-"-232 



_12 r 5 J T -±^Am - 9 A 2 r 3 + 40 A 2 r 2 m - 53 A 2 m 2 r 



+8mr 8 - 16 m 2 r 7 + 2AAm 2 r 4 J r _±2^L + i 8 A 2 m 3 I l/4r" 2 



r 3 + A< 



— r + 2 m 



— r+2 to 



~ 4 W 2Ar + 4Am 



-i 



(127) 



Averaging these expression in a spherical shell 2m < r < 3m, we have checked 
numerically that for A ~ 14(GM) 3 , 



(gy 

1-^232 / 



«1, 



(128) 
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clearly showing the possibility that in a shell of radius 0(m) around a "worm- 
hole" the "unstable" mode screens completely the classical curvature induced by 
the "wormhole" itself. In view of the general result that the zero-curvature met- 
ric realizes the minimum energy density of the classical gravitational field, the 
maximum screening of the curvature existing around of a "wormhole" allows the 
"unstable" mode of the quantized gravitational field to produce the minimization 
of the ADM (classical) energy that characterizes the Ground State of Quantum 
Gravity. 

7 The gas of wormholes: a possible ground 
state of quantum gravity 

Let us now, in this concluding Section, try to distill the main points and results 
of this work. In order to have a better understanding of what has been achieved, 
let us first describe again the strategy we have adopted to probe the stability of 
the classical ground state - the flat Minkowskian background - subject to the 
fluctuations of the quantized gravitational field. 

For definiteness' sake, for the classical solution of matterless gravity, our 
Ansatz for the classical background of QG, we have chosen the Schwarzschild 
solution but other possible starting points could be the Reissner-Nordstrom or 
the Kerr solution. From the classical standpoint such a background cannot be a 
good model for the vacuum, for it is well-known that the ADM0 mass of such 
solution is positive (Eadm = M), whereas Eadm = for flat space-time. But 
for the quantized gravitational field things may in principle be quite different, 
for Eadm receives contributions not only from the classical background field but 
also from the quantum fluctuations described by the quantum field related to 
the metric field by: 

9ij = r$ + h ij: (129) 

ryf^ being the Schwarzschild solution, fully characterized by the ADM-mass M 
(and Schwarzschild radius r s = 2GM = 2m). And it may happen, as it has been 
discovered in ref.0, that the full energy E(M) classical plus quantum, of the 
configuration (|129|) may turn out to be smaller than quantum energy E(0) of the 
Perturbative Ground State (PGS), where rjfj — > 5ij. If this indeed happened one 
would have finally got rid of the embarrassing PGS and of the unrenormalizable 
perturbative field theory that inevitably obtains upon it. 

As stressed in the Introduction, our hope that such strategy might lead to a 
similar situation, thus giving us a sensible QG, was based on a similar analysis 
performed on a non-abelian gauge theory, the SU(2) pure Yang-Mills theory (akin 
to QCD)0, where it was discovered that the energy E(B), of the quantized field 
fluctuating around a classical solution, describing a constant chromomagnetic 
field B, was in fact lower than the PGS-energy E(0). In the course of this analysis 
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it was also discovered that around this non-trivial classical background there 
exist "unstable" modes (whose frequencies of small oscillations are imaginary) 
that profoundly modify the structure of the theory even at very high frequencies, 



fatally undermining the paradigmatic Asymptotic Freedom ||12||. Of course, our 
hope was strengthened by the observation that QG may be looked at as a non- 
abelian gauge theory whose group is the Poincare group, acting upon the tangent 
Minkowskian spaces. 

The analysis reported in Sections 4 and 5 reveals that around an isolated 
"wormhole" : 

1. there exist an S-wave "unstable" mode [see (0), fl96|) and (|97j)1; 

2. the high energy "stable" modes are red-shifted with respect to the 
gravitons of flat Euclidean space, and realize the large energy gain 
AE(M) of eq.(El). 

As a result the PGS is clearly seen to be unstable: the quantized field finds 
it energetically advantageous to "concentrate" in a "wormhole" of radius 2GM 
and "mass" M, which redshifts the zero-point modes of the gravitational field, 
lowering in this way their energy density. 

But, as shown in Section 6, the single "wormhole" configuration certainly 
does not minimize the energy of the gravitational field, for the appearance of 
several other "wormholes" of mass M at an appropriate inter-wormholes (average) 
distance a - a gas or a lattice of "wormholes"- leads to a decrease of the energy 
density [see eq. ( |121| )1 provided, however, that a is bigger than the distance d = 
AGM, at which two of them coalesce into a single "wormhole". Taking into 
account such constraint, from ( |121| ) and ( |122|) we derive that an ensemble of 



"wormholes" with mass M ~ G~^ = m p at the (average) distance a ~ G^ = l p 
minimizes the part of the energy density that takes no account of the (negative) 
classical energy density stemming from the so far neglected "unstable" modes. 

At the end of section 6 by a numerical analysis we show that the variational 
amplitude A of the "unstable" mode [see eq. ( |124|) 1 can be chosen in such a way 
that the classical part of the Riemann tensor, averaged over a spherical shell 
(2m < r < 3m) around a "wormhole" is much smaller than the averaged Riemann 
tensor of the Schwarzschild solution. In this way the energy density can be further 
lowered, obtaining a state of the quantized gravitational field, fluctuating around 
a gas of wormholes of Planck mass with Planck distance separation, whose energy 
density is way below the energy density of the PGS: an excellent candidate for 
the GS of QG. 

Why do we believe that this latter statement is of significance for our un- 
derstanding of the fundamental laws of physics? First of all we must note, as 
emphasized in the Introduction, that the gas of "wormholes" that minimizes the 
energy density of the quantized gravitational field realizes the remarkable intu- 
ition of J. A. Wheeler that at the Planck distance continuous space may dissolves 
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in a kind of foam whose voids have the Planck size l p . Even though Wheeler's 
foam can only be seen when our space resolution reaches l p , just an impossible 
dream for presently available "microscopes" (the high energy accelerators), its 
effect on the structure and self-consistency of the QFT's of the Standard Model 
and for QG itself is enormous. In fact the space-structure upon which the observ- 
able quantum fields are defined turns out to be essentially discrete, comprising 
the Planck size space domains located at the interstices of the Planck size WH's 
which form the "points" of a "Planck Lattice" (PL). 

This latter statement might appear somewhat surprising, in view of the fact 
that the idea of a "foam" seems to retain a basic element of continuity of space- 
time which is lost in a lattice formulation. However, if we now look at the physics 
of the phenomenon of "WH-condensation" which, we have shown, leads possibly 
to the real GS of QG, we realize at once that no meaning can be attached to the 
value of the gravitational field or, for that matter, of any other physical field, in 
space-time domains whose (linear) size is less than the Planck length l p . And this 
for two reasons, first domains of such size are hidden beyond the horizen of the 
single WH and cannot affect the physics of the observable space-time domains (of 
the same size) apart from their gravitational effects, that are however compeltely 
screened by the unstable modes; and second the gas-like structure of the system 
of WH's with their fast quantum fluctuations cannot but "average" any "local" 
field over Planck-sized space-time domains, thus rendering unobservable, indeed 
unphysical, all field modes whose wave-lengths are smaller than l p . Without such 
modes in the spectrum of any physical field it obviously follows that there is no 
way to probe space-time deeper than l p , thus allowing us to effective describe 
space-time as an essentially discrete structure, a random lattice of average lattice 
constant a — l p . 

Even without going into the details of the significant advances in a fundamen- 
tal understanding of the SM that a research program based on the PL achievesQ, 
we can realize at once that the discretization of space in a PL introduces a fun- 
damental momentum cut-off A ~ m p , the Planck mass, thus removing at a basic 
level all divergences of QFT's, which since their formulations at the end of the 
20 's, have been a devastating conceptual difficulty, that the "conventionalistic" 
way out of "renormalization" could only throw "under the rug", arousing the 
anger of great minds such as Dirac's. With a PL, a faithful and adequate rep- 
resentation of the structure of the non-perturbative GS of QG, the ghosts of 
renormalization theory, such as the Landau's ghost, disappear beyond the hori- 
zon of the WH's, leaving us a decent perturbative expansion for the QFT's of 
the Standard Model, whose logarithmic divergences get converted to small cor- 
rections once one sets A = mjp. But this is not all, also the observable part 
of the quantized gravitational field, living on the PL, is cut -off at m p , and the 

7 A concise bibliography is to be found in ref.Jl3[ 

8 This obviously does not happen for scalar theories, such as the Higgs model, where there 
appear "quadratic" divergences in A 2 , which retain their devastating character, thus exposing 
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"nonrenormalizable" divergences 0[(GA 2 ) n ] become small corrections once the 
providential powers of 27T in the denominators, stemming from momentum inte- 
grations, are fully taken into account^. In this way, QG leaves the "limbo" of 
the theories of uncertain "self-consistency" to assume the status of a perfectly 
well defined QFT, with a well-defined perturbative S-matrix, whose quantum ef- 
fects are governed in a well-defined fashion by a very small coupling constant 
G, the Newton constant. And in view of all this one needs only the "Ockham's 
razor" to cut away much of the sophisticated theoretical developments, based on 
"super strings" , of the last 15 years. 

Let us end this paper by asking ourselves whether there remain open problems. 
Let us emphasize that we deem our results qualitatively robust, especially as for 
as the derivation of a quantum foam is concerned; an interesting open problem is 
then certainly the extension of the analysis to more complicated vacuum solutions, 
such as the Reissner-Nordstrom and Kerr metric. A more refined analysis, aimed 
at determining the detailed structure of the quantum foam and its approximation 
as a PL, is certainly needed. We hope to come back to this problem in the near 
future. Finally we think that "vexata quaestio" of the "cosmological constant" A c 
may have an answer within the framework of this paper, and that the "unstable" 
quantized modes of the gravitational field around each WH may have a lot to do 
with the surprising, and extremely welcome, negligible value that A c has in our 
Universe. 

We acknowledge the collaboration of S. Cacciatori and I. Spagnolatti in the 
early stage of this work. 



8 Appendix A: Christoffel Symbol and Ricci 
tensor for the background field 

We give detail calculations of some fundamental quantities generated by rj. We 
call 

/ Y^k \ 
r]= o" r 2 • (130) 

V r 2 sm 2 9/ 

For more convenience, we select the cartesiane coordinate: x 1 = r sin 9 cos (p, x 2 = 
r sin 9 sin <p, x 3 = r cos 9, and 

/ 1 2 3 \ 

/ 30 30 30 \ 

dx = — dr + x 3 cos (pd9 — x 2 dip\ — dr + x 3 sin ipd9 + x l d9; — dr — r sin 9d9 
\ r 1 r r 



their unplausibility. 

9 As we have been emphasizing in this work, the cut-off at m p fully justifies our neglect in 
the Hamiltonian of terms of order higher than two. 
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and then 



where 



As a result, 



ds 2 



v r J j r 2 

x ■ dx = ^2 x a dx a . 



r %a%b ( 1 
Vab ~ Oab + -JT \ ~ 2^ 



(dx) 



and, in particular, we obtain 



V r J 

We first calculate the Christoffel symbol, 

T 'ab S) = ^V rs (d a Vsb + d b rj sa - d s r] ab y 

Since 



2m 



x b x s 2m 



/ 2 \ X 

■ [3r — Amrj — 



and obtain, 



ab 



(s ) 2m S ab x r 



4m 



1 _ 2m /2m 



(^) 



2m 5 a fea; a mx a x b x r 



r 3 l _ 2m 5 (l_2m 
r y r 

In the coordinate (r, 0, </?) we have: 







r 1 - 

1 lm — 



r(r-2m) 

-(r-2m) 



3 + 








2m 



(i-H 







-(r - 2m) sin 2 9 J 



1 lm 



( o i 

1 r 

i 

r 

V o o 



o 
o 

sin 9 cos 9 J 



lm 



(00 \ 

cotgfl 
V \ cotg# 
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We can then calculate the Ricci tensor: 



dO) pr (s) _ -pr (s) p r ( s )pS (s) _ pr (s)p r ( s ) 
Jx afc 1 ra,6 1 a6,r ' L sa L rb L ab L rs 



Since 



we find: 



pr (s) 



m x a 

r 



pr 

ra,fe 



m5 ab ma;^ /g4m\ 

r / 



(139) 
(140) 

(141) 



and 



pr (s) 
a6,s 



2m 6m m 



3 + 



2m 



(l-^)J 



(5a 



■ t X\yX T 



\ fnx a Xf,x r x s 

~t~ O s lfX a X r ~t~ OsrXaXb J \ y 

From this equation, in particular we obtain 



15 + 



6r 5 — 10mr 4 



(142) 



m 



3 + 



2m 



(i_s)j 



+ 



mx a x b 



15 + 



6r — 10m 



2m 2 x a Xb 
3 (l-^) 



2 ' 



(143) 



and 



■pr (s)ps (s) 



2m 

5~0 sa^r 



mx a x s x r 



3 + 



2m 



2m mx r XbX s 
-^-d rb x b 



3 + 



r(l- 2 -f) 
2m 



2m 



m 



XaXb~ 



2m 



(l-^)J 



m 



2 ' 



(144) 



At the end, we have 



_pa (s)p S (s) 



2m mx a XbX r 

— OagXr 



2m 
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and result: 



R 



ah 



Owning to 



m x„ 
r 3 1 



2m 



m x a Xb 



r 



A 2m 



2m 2 



3 + 



>ab 



m5, 



ab 



4m 2 x a Xb 



2m 



3mx a Xb 



r 3 (l-^f) r 6(i_^) 2 r 5^_My 
2m 2 x a xi ) 

r 6 (l-^) 1 



+ 



x a x b m' 



r 6 1 



+ 



2m 2 5, 



2m\ 2 r 4 1 _ 2m 

—J 



+ 



m* x a x b 

r 6 ]_ 2m 
r 



2 + 



1 - 



2m 



m 

~ ab 



2m 2 x a Xf, 

r 6 ]_ 2m 



(I- 2 ?) 

x a Xb 2m 



m 3mx a Xb 

~Z^Vab + 



(l- 3 ?) 



2 i _ 2m ^ ab > 



the Ricci tensor of the background field is: 



u6 



//t / 3r 3r 

-7l — Vab + ~ Vab ~ ~ Oab 

\ Am. /m. 



r° v 2m 2m 

Going back to the coordinate (r, 9, if) in the matrix form, we have: 



R 



a (s) 



-7? 

25 



(145) 



(146) 



(147) 



In particular we find that the scalar of curvature is zero: R^ = 0. 



9 Appendix B: Decomposing the Regge- Wheeler 
matrix 

Using the technique due to T. Regge and J. A. Wheeler [||], we separate the angular 
part from radial part in the 3-dimension tensors <j>j(r, 9, <p), vector i^(r, 9, <p) and 
scalar r(r, 9, <p). 
The even part: 

<f> 1 1 (r,e,<p) = H(r)Y lm (9,<p), (148) 
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<j>j(r, 9, <p) = (d(r) + G 2 (r)d 2 e )Y lm (6, <p), 
<t>l(r, 9, <p) = ^Gi(r) + ^}dl + G 2 (r) cot 6d^ Y lm (0, <p), 
<fi(r,8,<p) = K(r)d e Y lm (e,<p), 



sm 

4>l(r, 0, <p) = < ¥£(de - cot 6)d v Y lm (6, <p), 
sm 9 v y 

^(r,0,p) = C/(r)y im (0,<p), 

v 2 (r,0,p) = V(r)0^m(M, 
U3 (r,0,<p)=y(r)d^m(0,<p), 

r(r,0,<p)=T(r)Y Jro (0,¥>)- 
#(r,0,<p) = O, 

02 ( r , e, v ) = F 2 (r)-^ (d e - cot e)d v Y lm (e, <p\ 

<l>l(r,0,ip) = -<ft(r,9,<p), 



The odd part: 



tf(r,9,p) = - I ^d v Y lm (9,(p), 
sm 9 



sm a 

<f>l(r, d^) = \^l {-^T-A + cot Ode - d*) Y lm (0, v), 
2 sine' Vsm 9 v J 

V!(r,9,(p) = 0, 

v 2 (r,9,p) = D(r) (—^) d e Y lm (9 , tp) , 
\ sin (7/ 

v 3 (r,6,<p) = D(r) sin9d e Y lm {9,tp). 
Obviously, the scalar has no odd part. 
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10 Appendix C: Explicit form of the eigenequa- 
tions of Hamiltonia 



We give the great details of calculations for obtaining the equations of Hamiltonia 
eigenfunctions and constrains. 
Laplacians of tensors 

We show how to operate the Laplaciano V 2 = ?? afe V a Vf, where the derivatives 
are taken covariantly with respect to the matrix 77^. 

For a tensor with rank 2, we find: 



V 2 /4 



j\ab - V 



^ab 



h 



j\a,b 



V 1 h r 
L rb lb j\a 



L jb IL r\a 



L ab'^jlr 



ab 



= T) 



h l A- ( h r T i — h i V r \ 4- Y i h r A-V i F r h s — F* V s h T 

n j,ab + \ a j L ra a r L ja ) b + 1 rb' L j,a + 1 rb L sa' l j 1 rb L ja'h 



— v r v 1 h s — V h 1 + r r v s h l —r r h i — r r r* h s 4- r r r s h l 

jb Sa r ih"'r n ~ L -ih 1 - ™"s 1 nM'ir nh 1 - xr lb 'i ~ nh L iV'".. 



= V 



ab 



n j,ab 



jb lb r,a 

9F* h r 
Zi ra u j,b 



jb ra s 



ab lb j,r 



Zi ja u r,b ^ u j L ra,b 



ab sr j 

h i T r 



ab jr lb s 



yi pr us 
1 rb L sa' L j 



_ r>pi -ps ir | -pr -ps u _rr n _ pr -pi is , -pr -ps u 
Zi rb L ja a s + 1 jb L ra a s 1 ab Ai j,r 1 ab L sr' L j + 1 ab 1 jr' l s 

For a covariant vector, we find: 



(167) 



V 2 B a = 7] bc B a \ bc = rf c IB, 



be 



>a\b,c 



^lc B r\b 



^bc B a\r 



be 



be 



( B a,b ~ Kb B s), c - Kc ( B r,b ~ F S rb B s ) ~ Tl c (B a , r - T S ar B s 

B a ,bc - Kb B s,c - Kb,c B s - Kc B r,b + r r ac r s rb B s - v r bc B a ^ r 



+ KeKM ■ (168) 

Using, in the place of T bc , the vectors calculated in Appendix A, We easily find: 



V 2 0} = 



2m 



r — 2m 



2 • 2 

r z sin 



+ cot 9d e - 4 



4(r-2m) 3 
3 ( t ) i + 2 



01- 



4(r — 2m) 



(<% + cot#)0 2 



r — 2m , 9 r — 2m ,o 

^0 2 + 2 —03; 



(169) 



V 2 2 = 



— % + h% + -y^-A 

r z sin T 



2r — 3m 



2m 



r 

2 cot 2 



d r + cot 



0^ + 



4(r — 2m) 



9 e 2 



4 nn ,q „r-2m ,, ^cot 2 6 1 
— cot ^^0^ + 2 — 0} + 2 — — < 



'I; 



(170) 
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r - 2m 2 1 2 1 
r r z r z sin 

r — 2m 2 cot 2 # 



-2- 



4 , + 4tr-2m)^3 



+ 7T- d^l + 4 cot 9 =— 0? 

r — 2m ,i cotq 2 9 , 9 
+ 2 — 0} + 2^-0 2 ; 



(171) 



V 2 2 = 



a (p + 2 f ^ 2 • 2 /) 



cottf 3r-9m cot 2 9 
+ ^^o g — 



- 2^,0? - ^ ^ _ 2 (mt , + ^2 + 2c^ 3. (172) 



V 2 0t = 



-4 



r r z r z sin 9 T 



r — 9m 



4r — 5m 3 cot 9 
a + — a 



2 



r 3 sin 2 6* 



a>} + 



2cot# 
r 2 sin 2 



r 3 sin 2 6> 



«9^3--(^ + 3cot^)03 ; 



(173) 



V 2 ^ 



3r - 4m 3 cot 2 0' 



! o2 2r — 3m 3 cot 6* 

a, h — s — a h — ^— a 



2(r - 2m) 



+ 



r 2 sin 6 

3 2(r-2m) 2 

^2 + 2-2/1 Py^i 

sin (/ 



o 2 cot 6^,0 2 cot „ , o 



r 2 sin 



r 2 sin 9 



(174) 



While for a vector: 

w = 



r — 2m 2 L 2 2r — m 2(r — 2m) 

a 9 I - 9 a 5 



Vl 



dg + COt 



v 2 - 



r 3 sin 2 



d v v 3 ; 



(175) 



V 2 t; 2 = 



^a - - + -a 



r r 
2(r - 2m) 



2cot# 



r — m 1 9 „ 
- — cot 2 e 



v 2 



r 2 sin 2 9 



d v v 3 ; 



(176) 
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V 2 T* 



-d 2 - -j + ^d r + - - - cot Ode 

ty ry Zi iy £ 



^3 



+ 



2(r - 2m) 



d,pVi + 



2cot# 



(177) 



Here with L we mean the usual operator of angular momentum expressed in 
polar coordinate. According to the decomposition of Regge- Wheeler given in the 
Appendix B, we obtain: 
The even part 

Substituting the even part in the Appendix B into eqs.( |169| ) -...( |177| ); taking 
into account the rules of commutation, we find: 



—L 



+ 



sin 9 
4 

sin 3 9 



(do — cot 9 



Y lm (e,(p) 



cot 9 
I 

sin 9 



-L 2 -2d 2 Y lm (6,(p) 



sin 9 



(d e -cote)Y lm {e,<p) 

cot 9 



Yl m (9,(p) 



d 2 e Y lm (9,<p) + - ir -Y lm (e,<p), 



sine "" v ,rj sin 6 9 
Y lm (0, f) = (l + cot 2 9)d e Y lm (9, ip) + 2 cos 9d[ 
- L 2 } sin 9de\Y lm (9,tp) = -2 cot 9 sin 91(1 + l)Y lm (9,<p). 



2/n;) v m ,> o,,,, az>2 Y lm(9, <p) 

v sin 3 9 



Defining, 



d, . jv* + - >!±V 



we easily obtain: 
V 2 $} 



A -4- 



r — 2m 



H(r) + A l(l + l)K(r) 



r — 2m 



V 2 $ 2 



D 



2G 1 (r)-G 2 (r)l(l + l) \Y lm (9,<p), 
r — 2m" 



Gi(r) + G 2 (r)«9 2 
r — 2m 



+ 2 ^if (r) + 4--^K(r)^ 



+2- 



Go(r) 



1(1 + l) + 2d 2 9 }Y lm (9,ip) 



(178) 
(179) 

(180) 
(181) 

(182) 



(183) 



(184) 



2/F>3 



r — 2m 



A -2- 



r — 2m . , 
+4 5— -RT(r) 



G 1 (r)-G 2 (r)(l(l + l) + d 2 
1 



sin 2 



2/1 d 2 + cot 



r — 2m . . 
+ 2 ^(r) 



--G 2 (r) 1(1 + 1) + 2d 2 \Y lm (9,<f), 



(185) 
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r — m n 2r — 9m 
A + 2^^<9 r — 



K( r ) + 1 if(r) - ^ d(r) 



+- 



/(/ + 1) + 1 G 2 (r) <9^ m (0,^) 



(186) 



V 2 $? 



r — m _ 2r — 9m 
A + 2^— <9 r 5 



/(/ + !) + ! 



sin^ 



#( r ) + ^ ^ _ 2 



(187) 



V 2 <^ 



A + 2 



r + 2m 



^(0+4^-^^)1-1 



sin 2 6 



(do 



- cot OWyYimtf,?), 



(188) 



r — 2m _ 9 /(/ + 1) 2r — m r — 2m 

■CC S 1 n «r - 2 5 



r 9 



t7(r) 



+ -l(l + l)V{r)\d e Y lm {e,ip) 



(189) 



r - 2m „, /(/ + 1) m n m 



r 

+2 r -^^B(r)}d e Y lm (e,v), 



^>2 y^«3 



C(r) 



(190) 



V 2 T* 



r - 2m a , /(/ + 1) m„ m 



r o 



^2 y^«3 



C(r) 



-2 =—5(0^^(5,^). 



(191) 



The Odd part 

Using eqs.( |178|) , (|179|) , ( |180| ), (|181|) and the expressions of the odd part in the 
Appendix B, we obtain: 

V 2 $} = 0, (192) 



V 2 $2 



sin^ 



„ r + 2m \ „ , s r — 2m „ . . 
A + 2 — F 2 r - 4 =— Fi r 



do — cot 6* 



d v Y lm (9,(p), 



(193) 



V 2 <^ 



„ r + 2m\ i r — 2m „ . , 
A + 2 — F 2 (r) - 4 —F 1 r 



sin 



<9 e - cot0 d v Yi m (6,ip), 



(194) 
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v 2 $ 2 



r — m „ 2r — 9m \ , . 
Z>, + 2^-0 r — J*i(0 



/(/ + 1) -2 



sin 6 1 



d v Yi m (9,ip) 



F 9 (r] 



(195) 



V 2 $? 



r — m „ 2r — 9m\ „ . „ 1 r, „ 
A + 2^^«9 r — Fx r - -3 W + 1) 



- 2 



sin 9 



(196) 



V 2 <^ 



„ r + 2m\ . „ r — 2m , . 
A + 2 — F 2 r - 4 =— Fx(r 



1 



2sin# 



sin 2 



2/3 V 



9! + cot 6d - dl 



Y lm (9,<p), 



VV = 0, 



r - 2m A , /(/ + 1) m „ m 



r 

r — 2m 



sin 9 



d^Y lm (9,ip) 



\sm0 ' 



(197) 

(198) 
(199) 

(200) 



11 Explicit form of eigenequations 

The eigenequations are given by: 



N 2 



8m 2 



r 3 (r — 2m) ^ 1 



+ A a + vj b + rfr = 4A< (201) 



where r is the scalar function explicitly expressed in eq.(|73l). Analogous to the 
decomposition of Regge-Wheller, we separate the radial part from the angular 
part. 



11.1 Equations for the even part: 

• Equations of eigen- vectors. 

Substituting the symbols of Christoffel and tensor of Ricci indicated in 
Appendix A into the explicit formulas of and v a in Appendix B for 
the even part, and using the formula calculated in C.l, we calculate: 
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the equation for H (a = b = 1): 

N 2 



r — 2m ^ 9 r — 2m „ 2 6r 2 — 27mr + 23m 2 
dt + 2 — & - - 



+ 4 



r 

r — 2m 



r 2 3 r 3 (r — 2m) 



H(r) 



+2- 



2m 



Z(Z + 1)^(0 
2G 1 (r)-/(/ + l)G 2 (r) \Y lm (0,<p) 



= i-H(r)AX+- 



r ~ 2m zTT( \ r ~3m V(r)l(l + 1) 
——d r U(r) - —^U{r) + — -— 



+iV : 



; Z(Z + 1) 



H{r)\Y lm {9,ip). 



(202) 



the equation for the combine function 2G ? 1 (r) — Z(Z + l)G 2 (r) (a = b = 2 
an then a = 6 = 3): 



iV 2 



r — 2m 



<9 2 + 2 



r — 2m 



<9 r 7 



(2Gi(r) 



, ,x 43r 2 - 12mr + 7m 2 

- id + DCW) + 5 r3(r _ 2m) gfr) 

-*—2—Kl + l)K(r)jY lm (9,<p) 
= {-4A 

+ iV 



2G 1 (r)-l(l + l)G 2 (r) 



2 Z(Z + l) r 



4 
3 



2d(r) - Z(Z + l)G 2 (r) 



17(r) 



2 r 2 



Z(Z + 1) 



the equation for G 2 (r)(a = b = 2 an then a = b = 2): 



(203) 



N 2 



r — 2m ^ 9 4r — 7m „ 2 

dl + — - — a + — 



G 2 (r) + 2 



2r — 3m 



2 ,9 2 + z(z + i)jr im (^v9) 

-4A2G 2 (r) + 2^ + iV 



,1(1 + 1) 



C 2 (r) 



2<9 2 + /(/ + !) ^ m (0 >¥ »). 



the equation for if(r) (6 = 3 e a = 1): 



(204) 



N 2 



r — 2m „ 9 4r — 7m n 
d 2 . + d, 



G 2 (r) + 2 2 -^—^K(r) 
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2d 2 +1(1 + 1) 



sin 2 9 



d v Yi m {9, if) 



V(r) 



A\2G 2 (r) + 2^ + N 



Go(r) 



2d 2 + 1(1 + 1) 



sin' 9 



d<pYi m (6,<p). 



(205) 



Finally for b = 3 and a = 2 we obtain the equation per G 2 . 



N 2 



r-2m 2 n r-2m _2_ 

°r Z „ 2 ar r 2 



+ 



4 3r 2 - 12mr + 7m 2 
3 r 3^ r _ 277^) 

-At- (do - cot 9)d v Y m (9, <p) 
sin y v 7 

|-4Af2Gi(r) - /(/ + l)G 2 (r)l + iV : 



2G 1 (r) - 1(1 + l)G 2 (r 

2r — 3m,„ _ . TJr , .) 
-—l(l + l)K(r)\ 



,1(1 + 1 



+ 7: 



l(l + l)G 2 (r) 

1K < r Vi) 



r — 2m 



r 



d r U(r) 



r — 3m 



r 



2G 1 (r) 
, U(r) 



1 



sin 2 9 



d e - cot 6)d v Y lm (6,<p). 



2 r u 

From this equation, we get the eq.(|79f) for the radial part. 
Equations of constrain are, 



(206) 



= and 
The second equation corresponds to 



TV 



0. 



m 



$v = — — — 

3\ l r ^ r _ 2m) 3 
The traceless condition leads to: 

2G 1 (r) - G 2 (r)l(l + 1) + H(r)]Y lm (9, (p) 
which coincides with flSOl). 



(207) 



For imposing the condition of transversality, we use commutating rules 
shown in paragraph C.l. 



Imposing the condition (|207|) for j = 1, one gets: 

Y lm (9,<p) = (208) 



d T H(l) + 3 . r l m , H(r) - K(r)l(l + 1) 



r(r — 2m) 
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- if j = 2 one gets: 

|r(r - 2m)d r K(r) + 4(r - 2m)K(r) + d(r) 

+ G 2 (r) [l - Z(Z + 1)] }0^ m (0, = (209) 

- Finally, if j = 3, one gets: 

|r(r - 2m)d r K(r) + 4(r - 2m)K(r) + <3i(r) 

+ G 2 (r) [l - Z(Z + 1)] }d v Y lm (e, <p) = 0, (210) 

and one verifies that the boundary conditions (|76|) and ( |75|) are equivalent to: 
{r 2 [#(r)£/(r) + Z(Z + l)#(r)V(r)]} | r=+00 



{r 2 [#(r)f/(r) + Z(Z + l)K{r)V{r)}} \ r=2m = 0, 
' (r) - O (J) 



r — > +cx) 



I V(r)-0(1). 

11.2 Equations for the odd part 

• Eigenequation: 



(211) 
(212) 



Putting 
Putting 



a 
a 



Putting 
Putting 



a 
a 



1, one gets identity 1 = 1. 

2, one gets the equation for F 2 (r): 

2 



N z 



A-3^ r +4W) 



+2- 



2r — 3m 



Fi(r) 



sin# 



d e - cot9 d v Y lm (0,<p) 



-4AF 2 (r) + —D(r) 



sin i 



cot# 



N' 



6 = 3 one gets equation (|213|) . 
1 and 6 = 2 one gets: 

r — 2m „ 2 r 2 — 5mr + 5m 2 



A + 2- 



^ — ^ 



r — 2m 



*i(r) 



(213) 
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+L(2-l(l + l))F 2 (r) 



sin 9 



d v Y lm (9, if) 



4AFi(r 
1 



sin 9 



Or 



d v Yi m (6, if) 



rj 



)D(r) 



(214) 



- Putting a = 1 and 6 = 2 one gets: 



N 2 



r — 2m „ 
A + 2 — <9 r 



2 r 2 — 5mr + 5m 2 



+- (2-1(1 + 1) )F 2 (r) 



r — 2m 
1 



F x (r) 



sin 



d e Yi m (8, (p) 



-4AF(r) - I (ft. - ?)l?(r) 



-<%y irn (6>,<^) 



-sin # 

Putting a = 2 and 6 = 3 one gets: 



(215) 



iV 2 



„ 2 \ „ , , 2r — 3m „ . . 
A - 3-<9 r + - F 2 (r) + 2 — F x r 



1 



2sin# 



dl + cot9d e -d 2 e 
sin 6* v 



-4AF 2 (r) + — £>(r) 



+ cot 9d e — d. 



2sin^ 



^d 2 
sin 2 6 f 



Y lm (0,f)- 



(216) 



These give (|55|), (B7|), which are used for the conditions of gauge. The 
traceless condition is automatically satisfied, while the (|207| ) gives: 



Per j — 1 identity = 0. 
Per j = 2: 



-d r F\(r)r(r — 2m) — 4(r — 2m)Fi(r) 



+F 2 (r) 



1(1 + 1) 



sin 9 



d v Y lm (8,<p) = 0. (217) 



Per j = 3: 



d r Fi(r)r(r — 2m) + 4(r — 2m)F 1 (r) 



-F»(r) 



1(1 + 1) 



—d v Y tm (6,<p) = 0. (21* 
sin 6 
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Finally the boundary conditions flT6| ) e ([75) become: 

l)r 2 F x {r)D{r) 



r=+oo 



r=2m 



and 



= "/(Z + l)r 2 F 1 (r)D(r) 
Di{r) — ► 0(1) per r — ► +00. 



(219) 
(220) 



12 Appendix D: Normalization of eigenfunc- 
tions 

According to the scalar product ([43]), one needs to normalize the quantity: 

r 2 sin 9 



Q 



0, <p)*\(r, 9, tp) + $f (r, 0, y?)$ 2 (r, 0, 



$f (r, 0, p)$|(r, 0, (p) + 2r(r - 2m)$f(r, 9, <^)$ 2 (r, 0, 
+ 2r(r - 2m) sin 2 0<&i*(r, 9, p)$l(r, 9, p) 
+ 2 sin 2 (r, 0, p)$ 3 2 (r } 9, <pj . 

We introduce: 

K = K(l); H = H(r); G x = G 1 (r); G 2 = G 2 (r); F x = F x {r) 

and 



(221) 



We indicate, 



where 



F 2 = F 2 (r);Y = Y lm (9,<p). 



Y lm {9,p) = —— JT(cos0)(-l)H- 



2ttC/ 



Zm 



P l m (cos9) = sin |m 
•Pi (cos 



2 (Z + m)! 



^ 2/ + 1 (/ - m) 



d(cos0)H 
1 d l 



Pi (cos 9) 



sin 2 ' 0, 



2 l l\dcos l 9 



with the conversion 



d0sin0Pz m (cos#)P; m (cos#) = C ( 



1 ni 
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and 



cLY = imY, 



d^Y* = -imY*. 



The even part is given: 



Q<+> = | 2 c/r^^^[$r^ + 2r(r-2m)($r^ + sin 2 ^f$?) 



+ 2sin 2 ^*$3 + $f 
r 2 sin # 



J drd6d/y^ H 2 Y*Y + 2r(r - 2m) (V 2 <9 e y*<%y 



^2 N 

+—^d in Y*d !n Y 



sin 2 6 v 
1 

'sin 2 9 

+ {G 1 + G 2 d 2 e )Y*{G l + G 2 d 2 e )Y 



+ -^1':.:. 2~ 



- cot 9)d v Y* (d e - cot 6)d v Y 



+ (d + C 2 -^9; + G 2 cot 00 fl ) F* 



+c 2 -^7i9 2 +c 2 cot^)y 

sin 6* ^ 



00 

dr 

2m 1 — 



-# 2 + 



2m 
r 

r 2 sin 



/• sin 

+ j^drdedip i _ ^ 



dr— 

2m 1 

2G 2 <% 



2m 



AT 2 Z(Z + l)2r(r - 2m) 



9 y* 

sin0 



sin# 



d^Y + 2G(Y*Y 



+ LYUYG 2 + 2GiY*G' 2 LF - 2G 2 <9^Y* 

-\-d% + cot Ode) Y 
sin # ^ / . 



(222) 



As a result, we have: 



/"CO ^ z f 

/ dr T7^ H 2 + ^ + l)^ 2 2r(r - 2m) 

J2m l\ z I 

+ 2G 2 + Z 2 (/ + 1) 2 G 2 2 - 2Z(Z + l)GiG 2 

+ 2G 2 2 / d0d<psm9 \d e (-^d^Y*) d e {^-£ V Y 



- <9 2 r* 



sin6> 



<9j + cot dde^j Y 



(223) 



We calculate: 



/ difdO sm6 I ^f^— d v Y] 
J v sm# ' 
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where 



And: 



m 

o 

m 



dxUl- x 2 ) [d x ((1 - x 2 ) m ^d m P l (x)) 



dx{l-x 2 ) -(m-l)x(l-x 2 ) m ^ 1 - 1 d m Pi(x) 



+ (l-x 2 ) — d m+1 P l (x) 

2 ,1 

dx 



m 



r 2 



m - 1) V(l - x 2 ) m - 2 d m Pi(x)d m P l (x) 

+ (l-x 2 rd m+1 Pi(x)<9 m+1 P,(z) 
- (to - - x 2 ) m ~ 1 d m P l (x)d m+1 P l (x) 

- (TO - l) 2 J( m ,l) + (TO - l) 2 if( m ,0 



+ — 2 (m — l)if( mji ) 



TO 



' J K0 = /\(1 - x 2 ) m - 1 d m P l {x)d m P l {x)dx 
K (m ,i) = J\(l - x^-^Ptix^P^dx 
H {m $ = Si x x{\ - x 2 ) m - 1 d m Pi{x)d m+1 Pi(x)dx 
- J\ x(l - s 2 ) ro 9 m P I (x)9 m+1 Pi(x)dx. 



i(m,l) 



dOdip sin # 



-<9 2 r 



2v * (-^^ 

V sin v 



dOdip sin 6* 



/(/ + 1)F* 



d?„ + cot 6d )Y 
1 



to j + cot #<%r 



sin 



1 



+ cot 2 6d s Y*d6Y - 2m' 



cotO 



+ 



sin 6* sin U 

m,Z) 



vrfK^X) + J(m+i,i) - m 2 J{ m ,i) - 2mH( m ^ + m 2 K( m>V) 

2to 3 j FG, 



C;,m+i 2 + m2 Ci m 2 — m 2 J( m n — 2ml 



+ 2m 2 H {ni)l) + 2m 3 J( mii) 



■ (m,Z) 



(rn,l) 



Summing (|224 ) and ( |226| ), ones K ets§ 

1(1 + 1) 



A 



(m,l) 



Cu 2 



+ 



r 2 



J m m(l - TO) + I m - Cl m 2 



2m 2 (m - l) 2 K m - Q M 2 + m 2 CiJ 



3 We omit the index I. 



(224) 



(225) 



(226) 
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+ J m m 2 (m — 1)(3 - m) + (1 + m 2 )J rn+1 

- 2m(m - l) 2 H m - 2m/J . (227) 

We calculate explicitly the integral ( |225| ) and obtain: 

I m = -I m + C lm 2 -2mj\l-x 2 ) m ~ 1 x 2 d m P l (x)d m P l (x)dx 

= -I rn + (2m + l)C lm 2 - 2mJ m (228) 

thus, 



J, 



±/. + til. 



m ± m 

m 



4~ — —Cl m 

m 



(229) 



And 



(1 - a; 2 )a m+1 P z (x) = 2mxd m P l {x) — (I + m)(l —m + l)d m ~ l Pi{x 



therefore 

-Zrrc 



-Cu 2 



J lm 



(I + m)(Z - m + 1) 



2m J-i 



2m 
1 

2m 

(Z + m)(Z — m + l)2m/ m _i 



(1 - x 2 )" 1 ^" 1 ^^)^ ^(^)^ 



Cj,m+i -(/ + m)(/-m + l)C / , 



(230) 



owing to Ci )in+ i 2 = Ci }m 2 (l + m + 1)(Z - m) e C^ m 2 = C^ m -i 2 {l + m)(Z - m + 1), 
one has: 

-Z m = G,m 2 + (I + m){l - m + l)I m _l 

= 2C ijm 2 + (Z + m)(Z + m - 1)(Z - m + 1)(Z - m + 2)I m _ 2 



2 2 (Z + m)! 

"1 Oz, m + 7: rrJo- 



(Z-m)! 



(231) 



As a result 



^ xP^dP^dx = -I + ^ [Pifdx - xPf{x) \\ 



1 



2Z+ 1 

(m - l)Ci )Tn 2 e J m = ^-Ci^ 2 . Now we have: 



(232) 



d m ~ l Pi(x)d (l-x 2 ) m ~'d m Pi 



dx 



2(m - 1) ^ d m - 1 P l {x)d m P L {x)x{\ - x 2 ) m - 2 dx 
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2\m-l 



- J* d m - 1 P l (x)& n+1 P l (x)(l - 
= 2(m - 1)^! - j 1 P(x)cr +1 p(z)(l - x 2 )" 1 - 1 
= 2(m - l)# m _i - [ y 1 m - 1 P,(a;)277ia;d m P,(a;)(l - 



2\m-2 



dx 



-(l + m)(l-m+l) J 9 m - 1 P,(x)9 m - 1 P,(x)(l-x 2 ) m_2 dx 
-2// m _! + (/ + m)(/ - m + 1) J m _i 



and then 



or 



(Z + m + !)(/ - m) _ J m+ i 
2 m 2 

(2Z + l)(Z + m + l)(Z-m) 
4m(m + 1) 



(233) 
(234) 
(235) 



Finally 



j m = J 1 (i-x 2 ) m+1 d m p l {x)d m p l {x)dx 



= -2 



X 



-1 



(1 - x 2 ) m - 1 d m+1 P l (x)d m P l (x)dx 



+ 2(m - 1) ^ x 2 (l - x 2 ) m - 2 d m P(x)d m P(:r) 
= -2P m - 2(m - 1) J m + 2(m - l)K m 



(236) 



which leads to: 



2m -1 Pf m 



2(m — 1) m — 1 



2/ + 1 



Therefore we find: 

A(,, ro) = Z(Z + 1 



4m(m — 1) 
2/ + 1 



2m - 1 + 



(/ + m + l)(/-m) 
m + 1 



(237) 



(m — 1) + m — / — m 



+ -m(m — l)(2m 



— 1)(2/ + 1) + 



1 m(m - 1)(2Z + 1)(Z + m + 1)(Z - m) 

2 m + 1 



- (Z + m + 1)(Z - m) + m 2 + -m(m - 1)(3 - m)(2Z + 1) 

2 



i 



2(m + 1) 
1 

2(m + 1) 



(2Z + 1)(Z + m + 1)(Z - m)(m 2 + 1) 



(m - 1) 2 (Z + m + 1)(Z - m)(2Z + 1) 



= --1(1 + 1). 



(238) 
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Thus one gets: 

Q« =/ — {H 2 + l(l + l)K 2 2r(r-2m) + 2G^ 
- 21(1 + l)GiG 2 + G 2 2 [l 2 (l + l) 2 - 1(1 + 1)] }dr 
Concerning the odd part, one gets: 



(239) 



Q { ~ ] = [ drd9d(psm9^— (2r(r - 2m) 
Jt, N 2 I 



1 

sinV 



F 1 2 ^d ip Y*d ip Y 



+F 1 2 d e Y*d e Y + 2F 2 2 -^- | (d e - cot 9)d v Y lm 



+ 1 2 F > 2 



sin 2 6 

l(l + l) + -^-dl + 2 cot Ode 
sin 9 ^ 



Y 



l 2 (l + lf 



= / dr—{2r(r-2m)F l 2 (r)l(l + l) + F 2 ' 

+ 2 ( m2 /^( l * ( ^ y)|2+i(, + 1) (-^ y ' r 

4 

+ cot 9d e Y*Y) (^j-Y*Y + cot 2 9d e Y*d e Y 
' v sin 9 



(240) 



Using the results from calculating the even part, we obtain: 



Q { ~ ] = 2r ( r - 2m W + WV) + F 2 



l 2 (l + lf 



-1(1 + 1) 



dr 

(241) 
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13 Appendix E: Study of eigenequations in the 
limit of flat space m = 

At the first we study the even part. The eigenequations and the equations of 
constrains becomes: 



d. + ~d r 

ry ry £ 



1(1 + 1) 4 



H(r) + -l(l + l)K(r) 



A 



+ - (2G 1 (r) - 1(1 + l)G 2 (r) = -A-^H(r 



1(1 + 1) 



d; + -d r - " ' ' + ., 

ry ry £ ry£ 



G 2 (r) + -K(r) = -4-^G 2 (r) 
r n 



dl + -a 



2 Q 1(1 + 1) 2 



r 2 



Gi(r) + -G 2 (r)Z(Z + l) 



(242) 



+ ^H(r) = -A^ Gl (r) 



dl + -d r 



1(1 + 1) 



ry ry £i ry 

+ (l - 1(1 + l))G 2 (r) 



K(r) + -H(r) - — d(r 



and 



2G 1 (r) -G 2 (r)l(l + l)H(r) = 



d r H (r) + 3 



H(r) 



K(r)l(l + 1) = 



(243) 



r 2 d r K(r) + ArK(r) + G 2 (r) \l - 1(1 + 1) 



0. 



We show that eigenequations are consistent with equations of constrains. Sub- 



stituting the first three equations of system ( |242| ) into the first equation of 

1(1 + 1) 



4A 



2G 1 (r)-l(l + l)G 2 (r)+H(r 
2G 1 (r)-l(l + l)G 2 (r)+H(r 



di + -d r - 



(244) 



Now we substitute the first and last equations of (|242|) into the second equation 
of ( |243| ) ; take into account the relations: 



3 r 



-d r H(r) 



-dr 



d r H(r) --d r H(r) 
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^d r H(r) 



d 2 H(r) 







-d r H(r) 



-H(r) 



d r H(r)\ - ^d r H(r) 
+ \d r H{r) - %H(r) 



-H(r) 



+ 3^(0 



and using the third equation of ( |243[ ) , we obtain: 

4A 

¥ 



d r H(r) + - /(/ + l)K(r) 



2 + -d r -1(1 + 1)-- 

d r H{r) + -H(r) - 1(1 + l)K(r) 
r 



(245) 



Finally, substituting the last three equations of system fl242f) into the last equation 
of (|243|) , we have 



r d r 



d 2 K(r) 



r 2 d r 



r 2 d, 
rd 2 K(r 

1 



d. 



r 2 d r K{r) - 2d r K(r) - Ard 2 K(r) 



-d r K(r) 


1 r 


r 2 d r K(r) 


r 


r 






1 








r 2 d r K(r) 



0: 



rK(r) - 2d r K(r) 
1 



~<9, 
r 



r 2 d r K(r) -2d r K(r) 



d r K(r) 



and the same for ( |243| ) , one gets: 
4A 



-3 r 



rK(r) K(r) 

J r 



r 2 d r K{r) + ArK(r) + G x (r) + G 2 (r) ( 1 - /(/ + 1) 







1(1 + 1) 



r 2 d r K(r) + ArK(r) 



+ G 1 (r) + G 2 (r)(l- 1(1 + 1)) 



(246) 



The equation ( [244j ) , ( |245|) and ( |246| ) show that the conditions of constrain ( |243|) 
are compatible with the eigenequation (|242|) . Thus, we can solve the eigenequa- 
tion by using the first equation of (|242| ) together with the eq.( |243|) . As a result, 
we have an equation for H(r) only: 



d 2 + 6 -0r+ 6 - l{ [ +1) 



H(r)=A±H(r). 



(247) 
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Introducing the variable p = 2^r and indicating " ' " as the derivative with 
respect to new variable, one gets: 



6 

H + -H' + 



1 (l + W-2) 



0. 



(248) 



SettingE 
one gets: 



H = p a J„(p), 



H' = ap a ~ 1 J u (p) + p a J' u (p) 
H" = a(a - l)p a - 2 J u (p) + 2ap a - 1 J' u (p) + p a J" v (p), 
and substituting it into (|248| ) : 



(249) 
(250) 



T/ . 2a + 6 
J "„(p) + J v + ( 1 



(J + 3)(Z-2)-6a-a(a- 1)~ 
? , 



Up) 



(251) 



with 2a + 6 = 1, so that the above equation is the equation of Bessel function 



with v 2 = (l + |) . Finally, 



From the second equation of fl243Q , one finds: 



(252) 



K(r) 



Ip 



-7/2 



2 ' + 3 



2y/\ 1 T/ 2\/Ar 

J 7-L.i 



ft r 5 /2 '+2 ft 

Recalling the relation for the Bessel function 

d 



(253) 



— ( Z -»J v ( Z ))=- Z -vj u+x ( Z ) ) 



we have 



Vx 



(254) 



11 </„ are the Bessel functions. 
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thus 



K(r) 



1(1 + 1) 



l + 1 J^(2^)- 2 -^'j,^(2^ 



-7/2 u l+i 



ft r 5/2 u l+i 



h 



(255) 



From the first equation of (|243|) one finds: 



2G 1 (r) - G 2 (r)l(l + 1) + cr^' 2 J l+ , U^f-) = 0. 



(256) 



From the last, instead, one find: 



Gx(r) + G 2 (r)[l - 1(1 + 1) 
-2. .ff (2^ r 



4c fl + 1. /„vT 

' 1(1 + 1) \ rs/2 J *+h r % 



hr 3 / 2 J l+ a 



h 



1(1 + 1) 



7 1 + 1 

2 r V 2 



4+a 2 



r 3/2 ft «+i /} 
(^ + 1) 2 



2ftr 3 /2 I fr 



r l/2 ^2 



/(/ + !) 

4 A 



,1/2 ^ 



2 



r 5/2 2+ 



ft 



**3/2 ^ 



ft 



ft 



At this point, we can adopt approximation ft — > 0, which is equivalent 
limit of J u+ i(z) for z — > +oo, 



(257) 
to the 



7T 



^-(^ + 1)2 



and one has: 



H(r) ~ (7 



\ 



tt V Ar 



(25£ 



AT(r) 



C 

W+i)\ 



ft /2 n VAr 

^Va 2 ^ cos 



4^ + 2 )f 



-2-n/A. 



i(Z + l) \ TT V A 



sm 



2— + 



(259) 
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2Gt(r) - G 2 (r)l(l + 1) 



h 2 l 



cos 



(1 + 1) 



7T 



(260) 



Gi(r) +G 2 (r) !-/(/ + !) 



ft 2 4A 



cos 



\/Ar 



/(/ + l)\l 7T V ArR 1 
From the last two equations, one finds: 



-(J + 3) 



7T 



(261) 



con 



G a (r) 



_ffi±iU(Z + l)\ 7T V A 



ft 2 



4A 



cos 



(262) 



It is clear then that the ( |101|) coincide with these solutions in the limit m — > if 
a(l) = i(J + l)f. 

Now we study the odd part. In this case, the equation becomes: 



6 



d; + -d r - 



1(1 + 1) 6 



dl + -d r 



1(1 + 1) 2 



+ 



= -4AF 2 (r) 



FAt) = -4AFi(r) 



F 2 (r) - -iq(r) 
r 



(263) 



Fo(r) 1 



l(l + lY 



(r 2 d r + 4r)Fi(r). 



(264) 



Again the conditions of constrain are compatible with the eigenequation. In fact 
one finds: 



-4A 



(r 2 d r + 4r)Fx 



(r 2 <9 r + 4r 
/(/ + !) 



a 2 + -a 



Z(Z + 1) 6 



2 

Z(Z + 1) 



# + -d r 

r 

4 

-Fi (r . 
r 



/(/ + !) 



F x (r) 
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Based on: 



r 2 a 



<9 r 2 Fi(r) 



r 2 d r 
r 2 <9, 



-d r F x (r) 
U(r) 



r^F^r) - 29 r Fx(r) - 4rd 2 Fi(r) 



-9 r [r 2 5 r F x (r)J -3d r F x (r) 



7 9, 



r 2 a r Fi(r) - -FAr) 
J r 



rdlFAr) = dl rFAr) - 2d r FAr 



-drFAr) 



-Or 



rFi(r) - -Fi 



one finds 



-4A 



(r 2 «9 r + 4r)Fx( 



'l - ^ *(r) 



a 2 + -a r 

r 



r 2 <9 r + 4r 



Z(Z + 1) 



F(r) - 



1 - 



1(1 + 1) 



Fo(r) 



(266) 



which show the compatibility of constrain. From the first equation of ( |263| ), one 

gets: 



As a consequence: 



1(1 + 1)' 



Fo(r) = c 



2V\ ( VaV 

+ r l/2^ J /+| ft 



5^_ / ^Arf 
2r 3 / 2 l+ * ^ ft / 



ft 



1 + 2 

r 3/2 °l+2 



V\r\ 2 V / Ar 1 / 2 



J^i I 2^- I - J z+ a 2 



ft 



ft 



In the limit ft — > 0, one finds: 
Fi(r) ~ c 



ft 12 1 / v/Ar 7/J n tt n 
^Va^ COS 2 — - !( ' + 1) 2, 



(267) 



(268) 



1 



i(i +iy 



Fo(r) 



-1- — 7-cos (Z + 2 - 

7rVA2 ft \ ft v '2 



ft /21 o VA . 
7r V A r ft 



vA r , . , . 7r \ 
2-^-/(1 + 1)5 . 



Again «(/) = — z(Z + 1). 



(269) 
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